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The weak-field expansion of the charged fermion propagator under a uniform magnetic field 
is studied. Starting from Schwinger's proper-time representation, we express the charged fermion 
propagator as an infinite series corresponding to different Landau levels. This infinite series is then 
reorganized according to the powers of the external field strength B. For illustration, we apply this 
expansion to 7 ^ i^i' and u ^ vy decays, which involve charged fermions in the internal loop. The 
• leading and subleading magnetic-field effects to the above processes are computed. 

o 
o 

>>: 
■ 

I. INTRODUCTION 

' Particle reactions taking place in the early universe or astrophysical environments are often affected by the back- 
ground magnetic field or the excitations in the medium . A typical example is the modification of neutrino index 
^ of refraction in the early universe or the supernova There one needs to compute the neutrino self-energy in the 
<^ \ medium or the background electromagnetic field or both. The neutrino index of refraction is then extracted from the 
■ modified dispersion relation of the neutrino. Another example is the plasmon decay 7* ^ z^z/ |l[ where the decaying 
photon acquires an effective mass through the effects of the medium or the background magnetic field. With such 
T-H : an effective mass, the above decay is kinematically permissible. Furthermore, the behaviors of electron propagators 
occurring in the internal loop of the above decay are also affected by the medium or the magnetic- field. This also 
leads to a modification to the plasmon decay amplitude. Finally, a more recent example is the enhancement of 
neutrino-photon scatterings due to the background magnetic field At the lowest order in the weak interaction, it 
-««^ is known that the amplitude for 77 — > vv is proportional to the neutrino mass . Hence the resulting scattering cross 
1-Q ' section is rather suppressed. On the other hand, the presence of the background magnetic field alters the structures of 
I internal electron propagators, such that 77 — > is non-vanishing at 0{Gf) even in the massless limit of neutrinos. 

Specifically, the 77 — > vv cross section is enhanced by a factor {mw /m)'^{B / Bc)^ due to a background magnetic field 
^ B d,^, where niw and m are the masses of W boson and electron respectively; Be = ni^ /e is the critical magnetic 
field. 

In the above processes, the relevant magnetic-field strengths are often smaller than the critical value Be- Therefore 
it is appropriate to expand the decay width, cross section or other physical quantities in powers of B / Be- In the 
literature, such an expansion is usually performed after the relevant amplitude is obtained For a more complicated 
Ci I process, it is not always convenient to do so since the amplitude to be expanded may be very cumbersome. In this 
article, we shall propose a more straightforward weak-field expansion, which is performed directly on the charged 
fermion propagator participating in the process. With the charged fermion propagators expanded, the physical 
amplitude can be easily expressed in powers of B / Be- To perform such an expansion on propagators, we shall begin 
with Schwinger's proper-time representation for a charged fermion propagator under a uniform background magnetic 
field M. It is useful to realize that Schwinger's representation can be recast into a series expansion in terms of Landau 
levels [pi. In the weak field limit B <^ Be, we shall demonstrate that one can reorganize the infinite series in powers 
of the field strength B. This is the expansion we are after. 

This article is organized as follows: In Sec. II, we will review Schwinger's derivation of charged fermion propagator 
in a homogeneous background magnetic field. Since the convention used by Schwinger differs from the currently 
popular convention, we shall repeat some relevant details of the derivation for clarification. We shall also illustrate 
how to rewrite Schwinger's result as an infinite series where each term is associated with specific Landau levels In 
the weak-field limit, we shall demonstrate how to rearrange the above series in powers of the magnetic- field strength B. 
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Finally, some technical issues relevant to the phase factor in Schwinger's proper-time representation will be discussed 
in this section. In Sec. Ill, we begin with a briefly discussion on our earlier work [Q, where the weak-field expansion 
technique is applied to 77 vD and its crossed processes in a background magnetic field To further illustrate 

the technique of weak-field expansion, we also calculate the decay rates of 7 — > z^i^ and the neutrino Cherenkov process 
J/ — » ^7 in a background magnetic field. Our results will be compared to previous calculations which are performed 
using exact charged- fermion propagators in the background magnetic field |10-lJ]. A few concluding remarks are 
presented in Sec. IV. 

II. CHARGED-FERMION PROPAGATOR IN A HOMOGENEOUS BACKGROUND MAGNETIC FIELD. 

A. The Exact Propagator Solution 

The Green's function G{x,x') of the Dirac field in the presence of a gauge field satisfies the following equation 

[i ^^e4-m)G{x,x')^5{x-x'), (1) 

where 5{x — x') is the Dirac's delta function and m stands for the mass of the Dirac field. We will follow the 
technique employed in Schwinger's paper which regards G{x,x') as the matrix element of an operator G, namely 
G{x, x') — {x' \G\x). Therefore, Eq. (m) may be written as 

(p-m)G-l, (2) 

with = Pf_i + denoting the conjugated momentum, which obeys the following commutation relations 

[n^,a;^] = ig^u, (3) 
[n^,n,] =ieF,,„ (4) 

with Ffj^i, = d^Ai, — dyA^ denoting the field-strength tensors of the gauge field. Eq. can be formally solved by 
writing 



(5) 



1 /"^ 
G=— = -i / ds(p + TO)cxp[-i(m^ - p2)s]. 

n-m Jo 

This integral representation for G implies that 

POO 

G{x,x'):^~i dse-'"''''{x'\{fi + m)U{s)\x) (6) 







where U{s) = e with H = -(p)^ = -11^ — ^eafj_^F'^'^ . We observe that U{s) can be viewed as the unitary 
time-evolution factor if one takes H as the effective Hamiltonian that evolves the state \x) according to 

\x{s))^U{s)\x{Q)), (7) 

where s is the proper time variable. One can now rewrite G{x, x') as 

G{x,x') = -i dse-'"'''[j^'{x'{0)\nf,{s)\x{s)}+m{x'{0)\x{s))], (8) 



where we have assumed n^(s) operates on \x{s)) and n^(0) operates on |a;(0)). We note that the operators and 
satisfy 

^ = -z[x„i/] = 2n„ 

^^ = -^[^^,i/] = -2eF^.^^ (9) 

for a constant field strength F^,^. In the matrix notation, we may write dx/ds — 211, and dYl/ds — —2eFTi. Further- 
more the transformation function {x' {Q)\x{s)) can be characterized by the following equations: 
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{idf, + eA^ix)) (x'mxis)) = (x'(0)|n^(s)|x(s)), 
-idl + eA^{x')) (x'mxis)) = (x'(0)|n^(0)|x(s)), 



(10) 



with the boundary condition: (a;'(0)|a;(s)) S^{x — x') as s 0. To evaluate Eq. (||), we first solve Eq. (H) and 
obtain 



n(.) 



x{s)-x{0) = (1 



n(o), 

-2eFs\ 



(11) 



This solution implies 



^ {x{s) - x{Q))K(x{s) - x^Q)), 

where K = i(eF)^ sinh"^ eFs. Therefore, one has 

1 i 
{x\^)\H\x[s)) = --^eaF - (x - x')K{x-x') - -tr(eFcotheFs). 

With this result, one can solve the first equation in (Eoh, which gives 



(12) 



(13) 



(a;'(0)|a;(s)) = C(x, 2;')s~^ exp 
X exp 



'-tr ln[(ei^s)"^ sinh(eFs)] 



^{x- x')eFcoth{eFs)ix ~ x') + ^ea^,^F'"'s 



(14) 



The factor C{x,x') can be determined by substituting Eq. ( p^ ) into the second and third equations in (|T^. Since 
the r.h.s. of these two equations are given by 



(a;'(0)|n(s)|a;(s)) = - [eFcoth(ei^s) - eF] {x - a;')(a;'(0)|x(s)), 
(x'(0)|n(0)|x(s)) = i [eFcoth(eFs) + eF] {x - a;')(a;'(0)|x(s)). 



(15) 



one then arrives at 



1 



idfj, + eAf^{x) - -eF^y{x' - x)" 



Ad'^ + eA^{x') + -eF^,{x' - xf 



1 

2' 



C{x,x') = 0, 
C{x,x') = 0. 



(16) 



Therefore C{x, x') is found to be 

C(x,x) = C'(x')exp 
= C(x) exp 



e£de (^A^ + ^F^A^-x') 
£de I^A, + ^F,^{^-xr 



(17) 



Here C'{x') and C{x) denote integration constants in x' and x respectively. Note that the integral A^ + ^Ff^^{(, — x'Y 
is a total derivative in the presence of a homogeneous field if the first homology group of the space-time M is trivial, 
i.e. Hi{M) = |l3|. Hence the phase factor is independent of the integration path connecting x and x' . One can 
further show that C{x') = C'{x). Therefore C(x') or C'{x) has to be a constant independent of x and x' . This 
constant can be determined by applying the boundary condition {x{s)\x' {Q)) — > (5'*(a; — x') as s — s- 0. One obtains 



C = -i{^TT)- 



(18) 
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with the help of the identity 



■ 2 2 

2" ^ dx = 



From Eqs. (|), (|T|)> & and (|l|), one arrives at 
where 



(19) 



(20) 



g{x,x') = -(4^)- 



ds 



m + -7 • (eFcoth(eFs) — eF) {x — x') 



X exp 



1 i 
--trln[(eFs)"^ sinh(eFs)] - -{x - x') (eFcoth(eFs) ) (a; - x') 



$(a;, x') EE exp |ie y d^'' 



M ^Ffj_y{^ X 



(21) 
(22) 



Note that the translation invariance is broken by the phase factor ^{x,x'). Note also that the phase factor $(x,x') 
vanishes if the path connecting x and x' is chosen to be a straight-line. In addition, if the background gauge field is 
a homogeneous magnetic field such that F12 — —F21 — B, one can show that 



a.yF^'' = 2Fi2a3 = 2Fr, 



1 B 

exp[ — trln(F"^ sinh F)] = , 

^ 2 ^ sinB' 

7(Fcothi^ ~F)x^ (7-a;)|| 



B 



sin B 



<J3 
as 



(7 • x)±e 



;-Fi2 0"3 



(23) 



a; {F coth F) x = a;| — S cot i? x']_ , 

with (a -5)11 — a°6° — and (a • &)_l = 0^6"'^ + a^b^ for arbitrary 4-vectors and 5^. Hence a| 
a']_ = a^a^ + a^a^. To simplify the notations, we shall denote (7 • p)\\{±) as 7 • P\\(±)- From the relations in (p3|), the 
propagator function C/(x,x'), which respects the translation invariance, becomes 



ds eBs 
sin(ei3s 



■ exp(— im^s + ieBsa^) 



X exp 



— (xm — eBscotleBs) x^^) 
4s " 



eBs 



sm{eBs) 



exp(— leBsfTa) 7 • xj 



(24) 



B. Weak Field Limit 



We find it is more convenient to cast (gj) in the form J 

d^p 



with 



g{p)= J d^xe'P^gia 



-ip(x — x ) 



Gip) 



(25) 



ds 



cos{eBs) 



exp 



2 2 , tan(eSs) 2 



IS I m — + 



exp(iei3s(T3)(m + 7 • pii ) 



eBs 

1-P± 
cos(e-Bs) 



(26) 
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One can further show that 

Sip) 



ds 



■ exp 



2 2 

-IS I m — P|| 



Q co8{eBs) 
cos{eBs) + 7i72 sin(e_Bs)] (m + 7 ■ P||) 

when the fohowing identities are apphed: 



tan(ei3s) , 



cos{eBs) 



exp(iZ(T3) = cos zl + i sin za^, 

(73 
C73 



<73 



«7i72- 



If we define a new variable v = eBs, then Eq. I27l) can be rewritten as 



1 



where 



g{p) = -i I dwexp(-zt)p)— [(rn + 7-p||)/i +7i72(to + 7-P||)/2 - (7-p±)/3] 



11 — exp(— iatanw), 

12 — exp(— iatanw) tanw, 
1 



I3 — exp(— iatanti)- 

with p = [m? — p^i^/eB and a = p\/eB. Because Ij{v) ~ Ij{v + mi) for j = 1, 2, 3, we get 

/ exp(— ii;p)/j = exp(— ipriTr) / dv exp{—ipv)Ij{v) 

Jo n=0 

1 



1 - e-v'^ 

It is sufficient to evaluate Ai since the other integrals are obtained using 

aa 

A:, = —{l-e-'n--Ai. 
a a 

To evaluate Ai = JJ^ (it;exp[~iatanw] exp(— ipw), we rewrite 



exp[— iatanw] = exp 



1 



(27) 

(28) 
(29) 

(30) 



(31) 



(32) 



(33) 



(34) 



The r.h.s. of this equation can be expanded using the Laguerre polynomials. Specifically, the Laguerre polynomials 
Ln{x) are generated by the following generating function 



(35) 



for \Z\ < 1. Upon multiplying Z on both sides of ( |35| ) and subtracting (|3^), one arrives at 

-xZ 



exp 



1 - Z 



= ^(L„(x)-L„_i(x))Z'^ 



(36) 



5 



where one sets = 0. Using the identity 

xZ + l 



cxp(-2j (37) 



I- Z 

with the identifications Z = — e^^*", x = 2a, and combining Eqs. (^2|), (|3j), and (^36|), one obtains 
Ai = / due"" ^(L„(2a) - i„_i(2a))exp(-2mi;)(-l)"exp(-ip'i;) 

= e""^C„(2a)(-l)" / duexp[-i(p + 2n)z;] 



n=0 



(-l)"C„(2a) 

n=0 

Using Eqs. (|30|), ( ^3|) and (|3^), one rewrites the propagator function G{p) into a simple form 



,e"(l-e-»ny ^'^ ■ (38) 



^ . f -^rfn(^)^ + <(«)i^ ^ ,7:|i^ (33^ 
where d„(a) = (— l)"e""C„(2a), d'^ = ddn/ da,p\ = ^'^'^ 

TO — Py 

D = (to + 7-P[|) + 7-p± 2 ' 

PI 

-D = 7i72(to + 7 -pii). (40) 

We note that, in the limit of extreme field strength, i.e. B ^ Be or B ^ Be, only part of the terms in Eq. ( |39| ) 
are relevant. In the strong field limit B ^ Be, only contributions from the lowest Landau level n — need to be 
kept. For the weak field limit B <^ B^, we shall demonstrate that the infinite series in Eq. ( |39| ) may be reorganized 
in powers of the magnetic field B. Therefore those terms with lower powers of B are more important in this limit. 
To reorganize the series, we first observe that 

-idnD + d'D I ^ -idnD + d'D 



E—lUn^J -t- _ ± \ ^ — '"r; 

pj + 2neB ~ pj ^ 1 + 2iifB 



Ti=0 n=0 



'2neB 



= E4MJ^ \-^DY,n%,ia) + DY,n'^dUa)). (41) 
^^qPl \ Pl / \ „=o 71=0 / 

The infinite series J2'^=o n'^ dn{ct) and J2'^=o ''^'^ '^ni'^) be evaluated with the the identity 

oo 

dnja) exp(— 2mf) = exp[— iatanw], (42) 

n=0 



which follows from Eqs. (34), (|3^) and (|37|). Let us proceed by taking a derivative d/dv on both sides of (|42|). This 
gives 



(— 2i)"'^ n"'^(i„(a) exp(— 2mw) = — exp[— iatanw] 



-la 

i"dn[a) exp[—2inv) = — 
cr 

n=0 

Taking this derivative k times, we find that 

oo ( / . \ 



(-2i)'=^n'=d„(a)exp(-2mi;) H — ^ + 0(a''-i) Uxp[-iatan?;]. (43) 



n=0 
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To be more specific, one can define U{v) = exp[— mtanw] following Eq. (p^. It can be shown that dyU = FU with 
F = —iaj cos^ V. Hence one can show that 



fc-i 



1=0 



[f'' + C^F'^-^dyF + C^F^-^dlF - 



ClClF^-\d„Ff] + K3(a) + Ki{a) + 0{a^-^) 



(44) 



Here = a\/\b\{a — fe)!] denotes the number of combinations of size b from a collection of size a. In addition, K3 and 
K4 denote the third and fourth derivative terms. They can be shown to be 



«;3(a) = ClF'^-^dlF + C^ClF'^-^d^Fd^F + C^C'^F^-^id^Ff 
K^{a) = C^F^-'^dlF + C^ClF'^-^d^FdlF + ClC^F^-^{dlFf 
+C^ClClF^-'^dlF{d^Ff + C^ClClF^~-^{d^Ff 



(45) 



Note that above formula for the expansion of d^U can either be proved by method of induction or can be read off 
directly from the combinatorial factor in the the expansion of [d^ + F)^ ■ 1 [Q. It is worthy pointing out that 
d^F(v = 0) = for all odd number k and the value of d^F{v — 0) when k is even can be computed directly. For 
example one can show that d^F{v = 0) = 2F{v — 0), and df;F{v = 0) = 16F{v = 0). Hence the order of a*^"^ and 



the order of a'^-^ terms read d^Uiv = 0) = 2Cf (-m)'^^^ ^ [^g^fc ^ 40C|](- 



Similarly, one can also show 



that the order of a'^ " term for the D term vanishes when n is an even integer while D term vanishes for all odd 
integer n. Hence, by setting = on both sides of Eq. (Esh, we obtain 



k /Q!\ 
2 V 2/ 



k-2 



n=0 
00 



(f) 



fe-4 



n=0 



fc-1 



k-2 



) + 2 








(f) 









k — 5 



(46) 



Here we only keep terms to the order of 0{a 



Since a = p^/eB, the leading terms on the r.h.s. of the above 



equation give up to order of 0{e^B^) contributions to G{p), as can be seen from Eq. (|4l|). Precisely we have 



n=0 



Pi 



2neB 



k=0 



-2eB 

pI . 



(f) 4^nf) 



D 



00 ^ 

00 

/£ = 



Pi 1 



w 



a 



fe-i 



4 

.2 X 



(f) 



k-i 



-iD 



i2C^D 



Pi 



eB 
D 



D 



Pi 

9 

Pl 



2 

Pl 

k-2 



fc-1 



iQ4ip) 



Pl 



eB 



2{k - 2)C^D 



eB 
Pl 



{pI? 



D 



(-ft) 



^eS + iQ2{p) + «^4(p) 



(p^ — m?)' 



;eB + it/2(p) + iGiip)- 



Pl 



fc-3' 



iQi{p) 



(47) 



Where (p) denotes terms of order e i? and e^B^ and iQi{p) denotes terms of order e B and e°B° . Therefore, 
by Eqs. (E3) and (p9), we arrive at 



iQijp) 



iD 



D 



p^ — m? [p^ — m?Y P 
i{j^ + m) 7172(7 -Pll + 'Ti) 



eB + i^^^ + ig2{p) + iQ^^p) 



p^ — rn^ 



eB + iQ2(p) + j^4(p)- 



(48) 
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The first term of Eq. ( fig ) is just the electron propagator in the vacuum, while the second term is its correction to 
0{eB). The corrections with higher powers in eB can be calculated in a similar way. For example, to evaluate the 
term Q2{p) and Giip): we need to compute a few identities. Note that C| — fc(fc — l)(fc — 2)/6 hence one can show that 
Sfe^o ^3 (~^)'^~^/^ = 1/(1 + x)'^ for all \x\ < 1 from consecutive differentiating the identity J2T=Qi~^)'' ~ ^/(^ ~^ ^) 



which is valid for all |x| < 1. Similarly, one can also show that X]fc^o('^ ~ ■^)^3 (^■^)'^ ~ 1/(1 + a;)^ — 4a;/(l 
all \x\ < 1. Therefore, one can extract the 0{e^B^) and 0{e^B^) terms from the series expansion in Eq. ( 
result leads to the following contribution: 



m. The 



«^2(W = T^hrD 



2e^B^ 



1 



4p] 



(p2 



(p2 



■,2)5 



where D and D have been defined in Eq. ([40|). Similarly, one can show that X^fc^o ^5 (~'^)'^~^ ~ 1/(1 + 
EZo C^i-^)'-' = 1/(1 + EZoik - ^)Cj{-xr-' = 1/(1 + xr - 6^/(1 + xy and j:Zoik - ^)Cl{-xr 
2/(1 + xY — 7x/{l + x)^ for all |a;| < 1. Hence the fourth and fifth order propagator iG4{p) can be shown to be 

-2plpl-3{piy] ^ , .,5ri5(pi)'-16pipi + 2pi 



(49) 



igi{p)^-[8ie^B^] 



D - Se^B" 



(p^ — m^y 



D. 



(50) 



C. Phase Factor 

In this subsection, we discuss how to treat the phase factor <i>(a;,x') as defined in Eq. (^). First, we note that 
$(x,x') is reduced to 



$(a;,x') = exp |ie J d^'^A^ 



(0 , (51) 



if the integration path connecting x and x' is a straight line. This choice of integration path is only for convenience 
since the integration in Eq. ( p^ ) is path independent provided that the vector potential Ap(^) is non-singular. Second, 
for a particular type of Coulomb gauge: 

A(e) = f (4-6,6 -a;'i,0), 

the exponent f^, d6^A'(C) vanishes, hence $(x, x') = 1. Therefore, by choosing the above Coulomb gauge, the phase 
factor ^(x,x') in the electron Green's function can be disregarded. Such a simplification is, however, no longer valid 
for more complicated processes where more than one electron propagators are involved in the process. To illustrate, 
let us consider an one-loop triangular diagram composed of three electron propagators. We denote vertices of the 
diagram as P, Q and R respectively. It is useful to recall that the full phase factor between two points P and Q is 



$(P, Q) = exp 




A^,{x) + ^F^,ix-Qr 



(52) 



according to equation (^2|). Here we use P^ to denote the coordinate of the point P. Similarly Q'^ and denote 
coordinates of the points Q and R respectively. As discussed before, one can set $(P, Q) = 1 by choosing the special 
gauge 

Aq(x) = Ao(x)+Aq(x) (53) 

with Ao(x) = B/2 ■ {—X2,xi,0) and Aq(x) — B/2 ■ (Q2,— Qi,0). Similarly, one can respectively set <l>(i?, P) and 
<f>((5, i?) to unity by choosing the gauges 

Ap(x) = Ao(x) + Ap(x) 

Afl(x) = Ao(x) + Ap(x), (54) 

with Ap(x) = B/2-{P2,-Pi,0) and Afl(x) = B/2 • (i?2, 0) respectively. Apparently, Aq{x), Ap{x), andylp(x) 
are distinct from one another. Hence they cannot be adopted simultaneously to set all phase factors to unity. In other 
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words, the phase factors shall give rise to a non-trivial contribution to the three-point amplitude. In fact this non- 
trivial contribution can be understood in an alternative view. Taking Eq. ( js^ ) as an example, the integrand on the 
r.h.s. of the equation can be written as A = A+ ^FAx where A = Afj,dx^ , A = A^dx^ , and FAx = F^^dx^{x — QY 
are all one- form. One can easily show that ,4 is a closed form, i.e., 

dA = 0. 



(55) 



Note that A is exact if the first homology group is trivial, namely, Hi{M) = 0. To be more specific, if the gauge 
function A^(a;) is regular everywhere, then the one- form A is also regular. Therefore there exists a zero- form oj such 
that A = dw \s an exact form. As a result, the line integration which defines ^{P,Q) is path independent. In our 
problem, we need to compute the product of three phases: <&(P, Q) ■ ^{R, P) ■ ^{Q, R)- It is then important to note 
that the one-form A in each of the above phases depends on the boundary point of the path, despite the fact that the 
gauge function A^{x) is regular. In other words, the gauge of A is chosen differently in each path, which then gives 
rise to a non-trivial phase for a three-point amplitude. Precisely one may isolate the boundary dependencies of A by 
writing, for example, A = A' — F^^dx^Q^ in the case of $(P, Q). Apparently, A' is an exact form universal to the 
three phases, while F^^^dx^^Q^ depends on the boundary point Q. Using this separation, one may rewrite each phase 
as 



$(x, x') = exp 



d£,^^(A^ + ]^F^,x'') 



exp 



-i-F 



d^,^'x"' 



(56) 



Let us denote the first factor ey-p[ier, dS^i'iA^ + ^F^^x")] as ^'{x,x'). Since $'(P,Q) • ^'{R,P) ■ ^'{Q,R) 
^'{Q,Q) = 1, we conclude from Eq. that 



$(F, Q) ■ <^>{R, P) ■ $(Q, R) = exp 



--{R-PYF^^{P-Qy 



(57) 



This is the nontrivial phase contribution one must attach to the amplitude of a three-point process when we write all 
weak field charged propagator according to Eq. (48). 



III. APPLICATIONS 



A. 77 vv 

The weak-field expansion derived in the last section has been applied to calculate the amplitudes of 77 — > vv and 
its crossed processes in a homogeneous magnetic field less than [Q. According to the discussion in the previous 
section, the magnetic-field dependencies of above amplitudes reside in two places: the first place is in the electron 
propagator which is affected by the external magnetic field, while the second place is in the overall phase which is a 
function of the field strength tensor F^^y. Let us now take 77 — > vv as an example for illustration. Since the incoming 
photon energies are much less than my/i we can calculate the scattering amplitudes using the following effective 
four-fermion interactions between leptons and neutrinos: 

G 

1^ = {vao,(X - l'o)vi) {ej^igv - gAl5)e) , (58) 

where gy = 1/2 + 2 sin^ 9^ and = 1/2 for I — e; gy — —1/2 + 2 sin^ 9^ and gA = —1/2 for / = ii,t. The Feynman 
diagram contributing to 77 — > vv is shown in Fig. 1 of Ref . [Q . We should remark that the contribution due to g^ 
is proportional to the neutrino mass in the limit of vanishing magnetic field. At 0{eB) in the limit B <C Be, it gives 
no contribution to the amplitude by the charge conjugation invariance. Therefore we shall neglect the contribution 
by gA- Likewise, we shall also neglect contributions by gv for I = ^,t, since —1/2-1-2 sin^ 9w = 0.04 <C 1. 

To 0(eB), the amplitude for 77 — > vv can be written as M = Mi -I-M2, where Mi arises from inserting the external 
magnetic field to electron propagators according to Eq. (|4|), whereas M2 comes from expanding the overall phase 
factor for the three-point function as shown in Eq. (p7|). Therefore one has 

yGpgv _ ^ ^ u [ dH 



Ml ^ iATiaeB^-^u{p2)i°'{l - 75)w(pi)ei e2 




(27r) 



7i72[7- (^ + fci)|| +-nA ^ i{y+m) ^ + m) 

[{l + kiY-m^f ^ P-'m?^ (/-fc2)2-m2 
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jii +m) ^ 7i72[7-^|| + m] ^ + ™) 



(/2-to2)2 

j^i +m) ^ i(J/+ m) ^ lil2[l ■ (^-^2)11 +H 

(/ + fcl)2-m2^ Z2_^2^ [{l^k2y 



(59) 



where gy — 1 — (1 — 4sin^ d^)/2 for and m is the raass of the electron. The first and second term in gy are the 
contributions from the W and Z exchanges, respectively. To write down the amplitude M2, we recall Eq. ( p7| ) which 
states that the overall phase factor for 77 —> vv is 



$(X, Y) ■ X) ■ $(y, Z) = exp 



-{Z-XrF^,{X~YY 



(60) 



With B in the forward z direction and choosing = (0, 0, 0, 0), we arrive at 

$(X,y) • $(Z,X) • $(y,Z) - exp |(^)(YiZ2 - >^2^i; 

ieB , 

^\-—{YxZ2-YiZx). 



(61) 



Since we calculate the amplitude only to 0{eB), the first term of the above expansion gives rise to Mi; while the 
second term gives rise to M2 which reads: 



M7 = i4:TTaeB 



Gpgv 



d^Y d^Z—{YiZ2-Y2Zi] 



dH d^q d^r 



12 ^^1^2 



X exp[-i(q - I - ki) ■ Y] exp[-i(r - q - k2) ■ Z]u{p2)^"{l - lb)v{pi) 

( , i(l/+m), ^iiej + m), ,i(/+m) 

X ^.(1 - 75),V^(--7.)^|T^(-^e7.)^5f^ 



P — rn 

We can recast the amplitude A/2 by using the equations 



(62) 



Y, exp[-i(g-Z-fci)-r] 



exp[-z(g -l-ki) ■ Y], 

ok 



Zi exp[— i(r — q — k2) ■ Z] = i— — exp[— j(r — q — k-)) ■ Z] 

an 



and the integration by part, such that 



M2 = iAnaeB 



Gpgv 
V2 



d^Y d^Z- 



i f dH d-^q d-^r 



(27r 



02 '^I'^a 



X exp[-i(q -l-ki)-Y] exp[-i(r - q - k2) ■ Z]u{p2)j°'{l ~ J5)v{pi) 
d d d d 
dli dr2 dl2 dri 



. J n ^ »(^+ m) . i(J + rn) . i{f + m) 
tr <j 7a(l - 75)1^ —i-'^^^t^)— Z:T(-*e7^) 



/2 — m2 



(72 — rn? 



(63) 



Before we proceed to compute Mi and M2, we wish to reiterate the validity of the above expansion. As we have 
pointed out in Ref. 0| that, by dimensional analysis, any given power of eB in the expansion of M is accompanied 
by an equal power of 1/to2 (for m > p) or 1/p^ (for p > m). Here p denotes the typical energy scale of external 
particles. Therefore, both eB/rn^ and eB/p^ are much smaller than unity for B Be = m^/e. Now performing the 
integration in Mi and M2, we obtain 



M=Mi+ M2 - ^^|^^S(p2)7a(l - 75)v{pi)J", 
v2v 47r 



(64) 



where [h5| 



J'^ = CiieiFe2)ik1 - k:^) 

+ C2[{eiFki){ki ■ e2)k^ + {e2Fk2){k2 ■ ei)fc?] 
+ C3[(eiFfci)e^ + {e2Fk2)e'i] 
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+ C4{eiFk2){h ■ e2)k1 + {e2Fkr){k2 ■ ei)k^] 

+ Cr,[{eiFk2)iki ■ £2)^2 + ie2Fk,){k2 ■ ei)k1] 

+ Ce[{e,Fk2)e^ + {e2Fki)e'^] 

+ Cr{k2 ■ • e2)[{Fkir + (^^2)"] 

+ C8(ei-e2)[(^^fci)" + (Ffc2)"] 

+ C9(fciFfc2)(ei-e2)(fc?-fc2") 

+ Cw{kiFk2){k2 ■ ei)(fci • e2)(fcr - k"^) 

+ Cii(fciFfc2)[(fc2 • ei)e^ - (fci • e2)e?] (65) 
with, for instance, {eiFt2) = £1^)^1/62 ^-^^d (Fki)" = F°'^kip. The coefficient functions C^'s are given as follows: 

Ci = — ^ f /[O, 0, 1] + /[O, 0, 2] - 4/[l, 1, 1] - 5/[l, 1, 2] + 2/[2, 1, 1] + 2/[2, 1, 2] 

+ t/[2, 1,2] + 2/[2, 2, 1] + 2/[2, 2, 2] - 5t/[3, 2, 2] + 2t/[4, 2, 2] + 2t/[4, 3,2]), 

C2 = — ^ f /[1, 1, 2] - 2/[2, 1, 2] - 3/[2, 2, 2] + 4/[3, 2, 2] + 2/[3, 3, 2] - 4/[4, 3, 2]) , 
V / 

C3 = — %:( 2/[0, 0, 2] - 4/[l, 1, 1] - 4/[l, 1, 2] - t/[l, 1, 2] + 2c/[2, 1, 2] + 2/[2, 2, 1] 

+ 2/[2, 2, 2] + 3t/[2, 2, 2] - 4t/[3, 2, 2] - 2i/[3, 3, 2] + 2i/[4, 3, 2] j , 
C4 = - ^ (5/[3, 2, 2] - 2/[4, 2, 2] - 4/[4, 3,2]), 

C5 = - A (^[1' 1' 2] + 2/[2, 1, 2] - /[2, 2, 2] - 10/[3, 2, 2] + 8/[4, 2, 2] + 4/[4, 3, 2]) , 

Ce = (2/[0, 0, 1] + 2/[0, 0, 2] - 4/[l, 1, 1] - 4/[l, 1, 2] - tl[l, 1, 2] 

- 4/[2, 1, 1] - 4/[2, 1, 2] - 2/[2, 2, 1] - 2/[2, 2, 2] + i/[2, 2, 2] 
+ 2i/[3, 2, 2] - 4i/[4, 2, 2] - 2t/[4, 3, 2]) , 

= ^ (/[1, 1, 2] - 2/[2, 1, 2] - /[2, 2, 2] + 4/[3, 2, 2] - 4/[4, 3, 2]) , 

Cg = 4^ f 2/[0, 0, 2] - 4/[l, 1, 1] - 4/[l, 1, 2] - i/[l, 1, 2] + 2t/[2, 1, 2] + 2/[2, 2, 1] 

+ 2/[2, 2, 2] + i/[2, 2, 2] - 4t/[3, 2, 2] + 2t/[4, 3, 2]) , 

^9 = - A (^[1' 1' 2] + 2/[2, 1, 2] + 4/[2, 1, 3] - /[2, 2, 2] - 10/[3, 2, 2] - 12/[3, 2, 3] 

+ 4/[4, 2, 2] + 4/[4, 2, 3] + 4t/[4, 2, 3] + 4/[4, 3, 2] + 4/[4, 3, 3] - 12t/[5, 3, 3] 

+ 4t/[6,3,3] + 4t/[6,4,3]), 

64 / \ 
Cio - ^ (/[4, 2, 3] - 4/[5, 3, 3] + 2/[6, 3, 3] + 2/[6, 4, 3] j , 

^11 = - A (^[1' 1' 2] + 2/[2, 1, 2] + 4/[2, 1, 3] - /[2, 2, 2] - 4/[3, 2, 3] - 4/[4, 2, 2] 

- 4/[4, 2, 3] - 4/[4, 3, 2] - 4/[4, 3, 3] + 4i/[5, 3, 3] - 4t/[6, 3, 3] - 4t/[6, 4, 3]) , (66) 

where 

/[a,5,c]= / / dy- — (67) 

with t = 2^j4a. Our result is an extension of the calculation in Ref. which considers only the low energy limit 
fci • A;2 <C rn?' . In such a limit, one can calculate M using the effective Lagrangian for 77 vv^ j9j and replacing one 
of the photons by the external magnetic field. 
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With the amphtude M, it is straightforward to compute the scattering cross section ctb (77 uv). in the back- 
ground magnetic field Since 77 — > vv could contribute to the energy-loss of a magnetized star, it is useful to compute 
the stellar energy- loss rate Q, which is related to as through [|6| 

1 [ 2d% f 2d% {ki-k2), , , , 

-(bJi + u}2)a-B[71 ^ i^i^)- (68) 



2(27r)6 J e'^IT - 1 7 e'^s/^ - 1 uji^2 

In Ref. Q|, Q is calculated based upon an approximated cross section obtained in the limit <C vn. Such a calculation 
is repeated in our earlier work |^ which is based upon the cross section ctb (77 vv) obtained from the amplitude M 
in Eq. (p^). We found that, for temperatures below 0.01 MeV, the effective-Lagrangian approach employed in Ref. 
I^works very well. On the other hand, this approach becomes rather inaccurate for temperatures greater than 1 MeV. 
At T — 0.1 MeV, our calculation gives an energy- loss rate almost two orders of magnitude greater than the result 
from the effective Lagrangian. Such a behavior can be understood from the energy dependence of the scattering 
cross section, as shown in Fig. 2 of Ref. It is clear that, for T — 0.1 MeV, Q must have received significant 
contributions from scatterings with lo k, m. At this energy, the full calculation gives a much larger scattering cross 
section than that given by the effective Lagrangian. By comparing the predictions of the full calculation and the 
effective-Lagrangian approach we conclude that the applicability of the latter to the energy- loss rate is quite 
restricted. While the effective Lagrangian works reasonably well with lo < 0.1m , it would give a poor approximation 
on Q unless T < 0.01m. 



B. 7 ^ uiy, u ^ wy 

In order to compare our approach with previous ones, we consider the simple two-body decay modes 7 — > i/i? [lO|-^T^] 
and V ^ lyj ||l^ in a background magnetic field. We shall limit the energies of incoming particles to be less than 
the pair-production threshold 2m. For if Ey,Ei, > 2m, the dominant decay modes should become 7 — + e^e^ and 
v — > i'e~^e~ respectively. For incoming energies below the pair-production threshold, it turns out that the photon 
momenta in both j viy and v ^ are space- like [|l8|. Hence the former process is kinematically forbidden. The 
amplitude of the latter process can be written as 

M{v{pi) ^ ^(^2)7(9)) = -^^e"u(p2)7^(l - l^Xpi) {gvn^p{-q) - 9Altp{-q)) , (69) 

where 11^^ and 11^^ are vector-vector and vector-axial vector two-point functions given by 

2 f d^k 



^lp{<D = -e" I [j^g{k - qh^^,g{k)] . (70) 

The factor Z is the wave-function renormalization constant of the photon field, induced by the effect of external 
magnetic fields. Since the deviation of Z from the unity is rather small, proportional to the fine structure constant 
a, we shall set Z = \ in our subsequent discussions. 

The structures of the two-point functions liap and 11^^ were given in previous literature P,p^ 

na/3(g) = A {qlg\\ap - 9itQ9||/3) + B {-q'ig^ap - q±aq±p) + C {q^gap - qaqp) , 
n^/3('Z) = C\\ (qlF^p + q\\c.{Fq)fi + q\\p{Fq)^^ 

+ Ci_ i-qlF^f) + q±c{Fq)f3 + q±|3{Fq)c^ . (71) 



We wish to remind the reader that = (q^)^ + [q^Y fo^' ^ magnetic field in the +z direction. The calculations of 
IIq^ and 11^^ for B < Be are straightforward using the weak field expansion derived in Eqs. (^8|), (|4^), and (|o|). Due 
to charge-conjugation and gauge invariances, the magnetic-field effects to IIq.^ and 11^^ begin at the order e^B^ and 
e^B^ respectively. The subleading contributions are then of the order e'^B'* and e^B^ respectively. The coefficient 
functions of Hap (q) are given by 
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A = 
B = 
C = 
+ ( 



7 , J5 ,2 . 26 52 w2 R ,4 

-45^^) +(315- 945;^^^^^)(^) + 



45 ^ 105 



1 • 2^^,S^2 



.45 
4 

105 



105to2--^%) 



44 
1575 



sin^ 9 



10 ^ 5 

2079 TO* '^Br' 



(72) 



where m is the photon energy while 9 is the angle between the the magnetic-field direction and the direction of photon 
propagation. For the coefficient functions of n^a(q), we find 



la 

'b^ 

ia 

'b^ 



70 TO^ 



■)^ + (- 



26 ^2 
945 to2 



r-sm"^ 6')( — y 

70to2 ^^B/ 



_8_ _86_^ . 2 

+ ^63" 945^'"' 



10 UJ^ . 4^, , B ^5 

jsm* 9)( — f 

693 TO* ^^Bc' 



10 ^ 4 B 5 

9){-) + 



693 TO* 



(73) 



It should be noted that the validity of weak- field expansion in Eqs.(^) and ( [73| ) also depends on the ratio r = 
up' SIT? 9 B^ / m? B^ , besides the requirement ^ 1- F^'' ^ sufficiently large photon energy such that r > 1, the 

expansion in Eqs. ( [72| ) and ( |7^ ) may break down. However, since we have limited the photon energy to < 2to, the 
ratio r is automatically smaller than 1. 

The computation oi v ~* width requires the knowledge of photon index of refraction n = —. The index of 
refraction can be calculated from the two-point function Hafiiq)- It is well known that n depends on the photon 
polarizations. For the magnetic field in the +z direction, the polarization states with distinct index of refraction are 
e'f^ = (0,0,1,0) and = (0, — cos6', 0, sinfl). Here we have adopted the convention that = (a;, cj sinfl, 0, cos0), 
i.e., photon propagates on the x — z plane with an angle 9 to the magnetic field direction. Hence e±_ is the polarization 
vector perpendicular to the x — z plane while eii lies on the x — z plane. The photon dispersion relation is given by 



ilia = 0, 



[l + iB sin^ 9 + iC] = iBiv^ siiP 9, 



[1 + iC] 



-iAlu^ sin^ 6', 



(74) 



where H^ = e^Hc/jef. Here the index a stands for the polarization states, namely a = _L or ||. Combining Eqs. ( [7^ ) 
and (uj), we arrive at 



(75) 



for polarization states a = _L and a =\\ respectively. Since the electomagnetic coupling constant is rather small, the 
left hand side of the above equations may be approximated by q^. Using the definition = uj'^ ■ {1 — rP) , we obtain 



«a = 1 + - 

TT 



ni, = 1 + - 

" TT 



1(J^)2 + (_ A + sin^ 9)i^^^ 

Ar^y^ ' inc ^ 135to2 ^ 



45' B/ ' 105 

7,5,2 , 13 26 ^2 2^, , B .4 

— ( — ) + ( \ ^sm^ 9)( — ^* 

90^Bc 315 945 to2 '^^ 



,B_ 

'Br. 



sm 



sm 



(76) 



It is seen that the leading contributions to n± and 7i|| agree with the results obtained by Adler |1^ . The next-to- leading 
contribution to n_L,|| depend on both the photon energy uj and the photon propagation direction |p^ . 

Given the above photon dispersion relation, we proceed to compute the v width in the sub-critical background 

magnetic field. We note that the most recent calculation oi v ^ width is performed by loannisian and Raffelt [pT| . 
Following their approach, we write the width of this process as 



r 



1 



167r£;2 



dcjY,\M\\ 

pols 



(77) 
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where A4 is the amplitude given by Eq. (|69|), Ei is the neutrino energy, and Wmax = niin(i?i, Wc) with Wc the critical 
photon energy beyond which the photon four momentum becomes time-like and the Cherenkov condition no longer 
holds. We note that, in deriving the above width, one has taken the collinear approximation that the particles in the 
initial and final states are all parallel with one another. The correction to such an approximation is small, proportional 
to the fine structure constant a. From Eq. (p9), we obtain 



\M\' = ^^e"(a)e-'(a) (pfpf + pf ^IpK'p. (78) 

where a =||,_L. The contribution by the vector- vector two-point function Ylap is negligible since both pi and p2 are 
approximately parallel to the photon momentum q and q^H-ap = due to the gauge invariance. The fact that both 
Pi and P2 are approximately parallel to q also has a consequence on the polarization dependencies of This is 

easily seen with 

'/'ns _ — 9 fr-n^^ ^ 2 



q"ltp = 2[C\\ql-C^qi){Fq)^. (79) 

For a B field in the +z direction, {Fq)a is nonvanishing only for a — 0,3. Given e'[ = (0,0,1,0) and eji" — 
(0, — cos6', 0, sin6') as stated earlier, one immediately see that \A4\^ vanishes for a photon in a _L mode. Hence the 
photon radiated from the neutrino is polarized, with its polarization vector lying on the surface spanned by q and B. 
The width of v ^ can be readily calculated using Eqs. (^7|), ( [7^ ) and ([7l|). We have 

r = ^aG^B^ ^.^g ^ r"""^ ^^Ei~cj)io^\C«-C^f. (80) 



Since C|| and C± are already given by Eq. ([73[), F can be easily determined once Wmax is specified. Since Ei < 2m, 
which implies cu < 2m, the photon refractive index is always greater than 1 as indicated by Eq. (|7^). Hence the critical 
energy tOc for photon dispersion relation to cross the light cone is greater than 2m. Thus LOmax = rnin(£^i,Wc) — Ei. 
The width F is given as follows: 



^^2G|^ 
135(27r)4 



1 S g / 8 1 Elsm^e\ B 8 



bO^Bc V105 49 m2 J ^ Be 



(81) 



Comparisons of our result with the earlier results of Ref. [ |12yi7|| are in order. First, we focus on the weak field region 
B < Be while Rcfs. |l^Jl^ considers the general magnetic field and the corresponding coefficient functions C'||_l are 
expressed in double integrals. Second, due to a different convention, the coefficient functions obtained in Ref. 
denoted as C'| ^, are related to ours via the relation |C[| — C± \ = 32^2^2 |C''| — 2C^| where 

,2 / j„ / J„.„~is(t>o n „.2^ 

(82) 

Jo J -I 
with 



Jo J-i 
C'j_^im^J ds J dve-'"''''R, 



9 1 — u 9 cos(eBsv) — cos(eBs) 9 
^° = ™ + -T'^^ + 2eBssHeBs) 

1 ~ V sm{eBsv) sm{eBs) — cos{eBs) cos{eBsv) , , 

R — 7r- . (80) 

si-a^ (eBs) 

To compare the two sets of results, it is useful to realize that one can rotate the integration contour, s —is, 
in the above integrals, provided qo = ui < 2m. In this way, the phase e~'**° turns into e"**" and becomes 
highly suppressed for a large s. For B < B^, such a behaviour permits one to simultaneously perform the 
weak-field and low-energy expansions with respect to C||.^. The results of expanding |C| — 2C^| may be orga- 
nized into the sum of the following series Y.n=o'^n{oj^ sm^ 9 /rn^YiB^ I BlY , Y.n=o^n{uj^ srn^ 9 / m^Y {B^ / BlY+^ , 
X]n=o Cn(w^ sin^ 6 / m?)"" {B^ / B^)""^^ ■ ■ ■. One observes that the coefficients ai and 02 correspond to the 0{uP'B^) and 
0(w^i3^) terms in our \C\\ — C\_\ respectively. We found that all the coefficients a!^s vanish. This is indeed reflected 
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in our calculations where the OiuPB"^) and 0{uj'^B'^) terms in \C\\ — C±\ vanish as well. We also found agreements 
between the coefficients 60,1 and the corresponding 0{B^), 0{uj'^B^) terms in |C|| —C±\. Although we did not compare 
the coefficient cq with the 0{B'^) term in Cy, due to the growing complexity in computing the general coefficients c-s, 
the agreements we just found with respect to the first two series seems rather compelling. Due to these agreements, 
we also confirm the statement made in Ref. [0 that the earlier calculation on 11^^^ is incorrect. 

From the above comparisons, we have seen that our approach, in spite of less general, is convenient for obtaining 
the analytic amplitudes of physical processes in a sub-critical background magnetic field. In such a magnetic field, it 
suffices to know the leading and sub-leading terms in the weak-field expansion. Our approach produces those terms 
directly from Feynman diagrams. 

The work on extending the present analysis to the more complicated processes, such as the photon splitting 7 — > 77 
and the pair production 7 —^ e^e~ is currently being pursued. For the latter process, we have exploited the analytical 
properties of the vacuum polarization function If^" in the background magnetic field. For a sub-critical magnetic 
field, it is possible to obtain a simple expression for the absorption coefficient(the pair-production width) for arbitrary 
photon energies . This is an improvement to the previous work where a simplified expression is possible only for 
Lu ^ m [ pl| . For the former process, 7 — > 77, our result shall serve as an additional check to the previous results 

IV. CONCLUSION 

In this paper, we have developed the weak-field expansion technique for processes occurring in a background 
magnetic field. This expansion is performed with respect to internal electron propagators which are affected by the 
background magnetic field. In some processes, our approach is valid for general external momenta even if they are much 
greater than the electron mass m. For external momenta much greater than m, the effective-Lagrangian approach is 
no longer appropriate. To illustrate this point, we calculated the amplitude of 77 — s- under a background magnetic 
field, and consequently determined the stellar energy-loss rate Q due to this process. It is interesting to find that the 
effective-Lagrangian approach is inappropriate for computing the stellar energy-loss rate due to 77 —^ vD, unless the 
star temperature is less than 0.01 m. This result reflects clearly the importance of our approach. In fact, our approach 
can be applied to many other processes. In this regard, we also discussed the processes j ^ i/v and 1/ ~* i/j under 
a strong background magnetic field. We found that the validity of weak-field expansion with respect to the above 
processes are also determined by the parameter r = uj^ sin^ 6B^ /rn^B^, besides the requirement B < Be- For energy 
below pair production threshold, the parameter r is less than 1, which causes no trouble to the weak-field expansion. 
We found that j vD is kinematically forbidden while 1/ —> 1^7 is permitted by the phase space. Our predictions 
on the latter process agrees with previous works ||l^. It has also been pointed out that our approach, although less 
general, is convenient for obtaining the analytic amplitudes of physical processes in a sub-critical background magnetic 
field. 

We are currently extending the weak-field expansion technique to the photon splitting process 7 ^ 77 |l^,|2^ and 
the pair production process 7 — > e^e~. ||l^ , pT[| . Both processes are of great interests in the physics of pulsars on which 
the background magnetic fields are close to the critical value Be- 
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